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Abstract We present the results of a statistical study 
of the star formation rates (SFR) derived from the 
Galaxy Evolution Explorer (GALEX) observations in 
the ultraviolet continuum and in the Ha emission line 
for a sample of about 800 luminous compact galaxies 
(LCGs). Galaxies in this sample have a compact struc- 
ture and include one or several regions of active star 
formation. Global galaxy characteristics (metallicity, 
luminosity, stellar mass) are intermediate between ones 
of the nearby blue compact dwarf (BCD) galaxies and 
Lyman-break galaxies (LBGs) at high redshifts z > 2 - 
3. SFRs were corrected for interstellar extinction which 
was derived from the optical Sloan Digital Sky Survey 
(SDSS) spectra. We find that SFRs derived from the 
galaxy luminosities in the far ultraviolet (FUV) and 
near ultraviolet (NUV) ranges vary in a wide range from 
0.18 Mq yr~ x to 113 Mq yr" 1 with median values of 
3.8 Mq yr _1 and 5.2 Mq yr" 1 , respectively. Simple re- 
gression relations are found for luminosities L(Ha) and 
L(UV) as functions of the mass of the young stellar 
population, the starburst age, and the galaxy metal- 
licity. We consider the evolution of L(Ha), i(FUV) 
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and L(NUV) with a starburst age and introduce new 
characteristics of star formation, namely the initial Ha, 
FUV and NUV luminosities at zero starburst age. 

Keywords Galaxies: irregular — Galaxies: luminos- 
ity function, mass function — Galaxies: starburst - 
Galaxies: star formation — Galaxies: statistics 



1 Introduction 



Cardamone et all |2009h first draw attention to galaxies 
at redshifts z = 0.112 - 0.360 which were named "green 
peas" because of their compact structure and green 
colour on the gri composite Sloan Digital Sky Survey 
(SDSS) images. Specific colours of these galaxies are 
mainly caused by the very strong [O in] A5007A op- 
tical emission line. The equivalent widths EW(A5007) 
of this line in "green peas" redshifted into the SDSS r 
band can be as high as ^lOOOA, resulting in a green 
colour on SDSS images. ICardamone et al.l (|2009l ) stud- 
ied a sample of 251 colour-selected galaxies. Some of 
the galaxies from this sample are active galactic nu- 
clei (AGN). However, most of "green pea" galaxies are 
found to be strongly star-forming ones with high star 
formation rates (SFR) of ~ 10 M© yr -1 . These galax- 
ies are characterised by low metallicity, stellar mass of 
M* ~ 10 8 - 5 - 10 10 M e , high specific SFR (SSFR) (up 
to ~ 10~ 8 yr -1 ) which place them between nearby blue 
compact dwarf (BCD) galaxies and high-redshift (z > 
2 — 3) UVTumin ous Lyman-break galaxies (LBGs, see 
Giavalisco 20021 for a review). The available Hubble 
Space Telescope (HST) high-angular resolution images 
of a few "green peas" reveal complex morphology on 
small spatial scales with several regions of star forma- 
tion and an exte nded stellar component likely consist- 
ing o f o lder stars ( Cardamone et al1l2009t Amorm et al 



20121 ). ICardamone et al.l |2009lU uggested that "g 
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pea" galaxies may be occurrences of the star formation 
mode prevailing in the early Universe. This galaxy class 
therefore may provide an excellent opportunity to un- 
derstand in great detail many processes under physical 
conditions approaching to those in high-redshift galax- 
ies. 

The oxygen and nitrogen chemical a bundances in 
star-fo rming "green peas" were studied by Amorm et al.^ 
( 2010). These galaxies are revealed to be genuine metal- 
poor galaxies with mean oxygen abundances of ~ 20% 
solar. The N/O ratios are found to be unusually high 
for galaxies of the same metallicity. Detailed study 
lead lAmorm et all |2010h to the conclusion that known 
general properties of "green peas" , namely high SSFR, 
extreme compactness and stellar mass, seem to be un- 
common in the nearby universe, suggesting a short and 
extreme phase of their evolution. The possible action 
of both recent and massive interaction-induced inflow 
of gas, as well as selective metal-rich gas loss driven by 
supernova winds are discussed here as main drivers of 
the starburst activity in "green peas" and their oxygen 
and nitrogen abundances. 

The first direct radio detection with low frequency 
Giant Metrewave Radio Telescope (GRMT) obser- 
vations and discussion of the "green peas" proper- 
ties comprising propertie s of a new class o f sub- mJy 
sources were reported by IChakraborti et al.l (|2012f ) . It 
was shown that this detection may imply large mag- 
netic fields ( > 30 //G) in "green peas" under reason- 
able assumption about cosmic ray diffusi o n and total 
energy consideration. IChakraborti et aL ( 20121 ) con- 
cluded that seed fields were amplified significantly (up 
to fiG) because of turbulence as protogalactic and sim- 
ilar structures formed. 

Detailed examination of a large sample of 803 star- 
forming luminous compact galaxies (LCGs) in the 
redshift ran g e z = 0.02 - 0.63 was carried out by 
Izotov et al.l (|201lf ). These gala xies were selected from 
the SDSS Data Release 7 (DR7) jAbazaiian et al.ll2009t ) 
and comprise a complete spectroscopic SDSS sample of 
strongly star-forming LCGs with reliably derived chem- 
ical abundances. Their global properties are similar to 
those of the star-forming "green pea" galaxies. How- 
ever, in contrast to "green pea" galaxies, the LCGs are 
selected on the base of the both their spectroscopic and 
photometric properties. Applied selection results in a 
~ 10 times larger sample, with galaxies spanning a red- 
shift range ab out ~2 times larger as co mpared to "green 
pea" sample (jCardamone et al ] l200fll) . For LCGs, the 
oxygen abundances 12 + log O/H are found to be in 
the range 7.6 - 8.4 wit h the median value o f ~ 8.11 con- 
firming the results by I Amorm et al. I fl2010h for a subset 
of the "green pea" sample of ICardamone et al] ([20091 ) . 



The ranges of oxygen abundances and heavy element 
abundance ratios in LCGs are similar to those of nearby 
low-metallicity BCDs. In the [O m1A5007/H/3 vs 



[Nn]A6583/Ha diagnostic diagram (jKauffmann et al 



2003) the LCGs are shown to occupy the region of 
high-excitation star-forming galaxies. The SFRs, de- 
rived from the Ha line emission in the LCGs are re- 
vealed to vary in the large range of 0.7 - 60 M© yr _1 , 
with a median value of ~ 4 M© yr _1 which is about 3 
times lower as comp ared to star-forming LBGs at z^3 
(jPettini et alj|200lh . The SSFR in LCGs is extremely 
high and it varies in the range ~ 10~ 7 - 10~ 9 yr _1 . All 
these properties imply that LCGs are likely the closest 
local counterparts of the high-redshift LBGs and Lya- 
e mitting galaxie s. 



Guseva et al.l (|2011l ) carried out the spectroscopic 
analysis of HG 031203, the host galaxy of a long- 
duration gamma-ray burst (GRB). The galaxy prop- 
erties such as the oxygen abundance 12+ log O/H = 
8.20 ± 0.03, extinction-corrected Ha luminosity L(Ha) 
= 7.27xl0 41 erg s"\ stellar mass M* = 2.5xl0 8 M Q , 
SFR(Ha) = 5.74 M yr" 1 and SSFR(Ha) = 2.3 xlO" 8 
yr _1 in HG 031203 are found to be in the range cov- 
ered by the LCGs properties. This fact implies that the 
LCGs with extreme star-formation, that also comprise 
"g reen peas" as a subcla ss, may harbor GRB. 



Pilvugin et al.l ([20121 ) analysed the oxygen and ni- 



trogen abundances derived from global emission-line 
SDSS spectra of galaxies using the direct method based 
on the electron temperature determination from the 
[O in] A4363/(A4959 + A5007) emission-line flux ra- 
tio and the two strong line O/N and N/S calibrations. 
Three samples of objects were comp ared, including the 
sampl e of "green pea" galaxies by ICardamone et al ' 



(200 9) with the dete cted [O m] A4363A auroral line. 



Pilvugin et al.l (|2012f) concluded that the high nitrogen- 
to-oxygen abundance ratios derived in some "green 
pea" galaxies may be due to the fact that their SDSS 
spectra are the ones of composite nebulae made up of 
several components with different physical properties. 

The local analogues of the strong Halpha Emitters 
(HAEs) dominated the z ~ 4 L BG population are iden- 
tified bv lShim fc Chary I (|2012h . Using the SDSS spec- 
tra authors show that at z < 0.4 only 0.04% of galaxies 
are classified as HAEs with equivalent widths EW(Ha) 
of > 500 A, comparable to that of z ~ 4 HAEs. Local 
HAEs have lower stellar masses and lower UV lumi- 
nosities than the z ~ 4 HAEs. On the other hand, 
their Ha-to-UV luminosity ratios and SSFRs are con- 
sistent with those of z ~ 4 HAEs indicating that the 
local analogues are the scaled-down versions of high- 
z star-forming galaxies. Compared to the previously 
studied local Lyman-break analogs (LBAs) of the z ~ 2 
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LBGs which were sel ected using rest-frame UV fluxes 
(|Heckman et al . 2005), the local HAEs show similar UV 
luminosity surface densities, but lower metallicities and 
lower stellar masses. This supports the idea that lo- 
cal HAEs are less evolved galaxies than the traditional 
LBAs. Local HAEs show a strong He II A4686 A emis- 
sion line in the stacked spectrum, implying a population 
of hot young (< 10 Myr) massive stars, similar to that 
seen in some Wolf-Rayet galaxies. The local HAEs also 
have properties similar to those of "green pea" galaxies. 

In present paper, we extend the study of the proper- 
ties of the "green peas" by fur ther analysis of about 800 
LCGs bv llzotov et all (1201 lh . selected from the SDSS 
DR7. The selection criteria of galaxy sample arc briefly 
described in Section [2] The correction of LCG fluxes 
for extinction is discussed in Section [3] In Section 2] 
we carry out the statistical investigation of dependence 
of galaxy luminosities on other LCG characteristics. In 
Section[5]we discuss star formation rates of LCGs. The 
SFRs are derived from the extinction-corrected lumi- 
nosities L(Ha), L(FUV) and L(NUV) in the Ha emis- 
sion line, the far ultraviolet (FUV) and the near ultra- 
bviolet (NUV) ranges from Galaxy Evolution Explorer 
(GALEX) observations. The luminosity function for 
LCGs is discussed in Section ??. We summarise our 
results in Section[B] We assume Ho=75 km s _1 Mpc -1 
for distance estimates. 



2 Sample selection, observational data 

It is noted in S ectionffl that LCGs, in c ontrast to "green 
pea" galaxies ([Cardamone et al.ll2009l ). are selected on 
the base of the both their spectroscopic and photo- 
metric properties. Selection criteria and LCGs sam- 



ple pr operties are described in detail by llzotov et al 



(|201lh . Briefly, these criteria are as follows : 

- the extinction corrected luminosity of the H/3 emis- 
sion line is greater than L(H(3) = 3xl0 40 erg s _1 ; 

- the equivalent width of the H/3 emission line is 
high, EW(H/3) > 50A. This criterion leads to selection 
only objects with strong emission lines in their spectra 
and thus the ones containing young starbursts with ages 
3 - 5 Myr; 

- only galaxies with well-detected [O in] A4363 A 
emission line in their spectra, with a flux error less than 
50 %, are selected. This criterion allows an accurate 
abundance determination using the direct method; 

- only the star-forming galaxies were selected. 
Galaxies with obvious evidence of Seyfert 2 features 
are excluded; 

- galaxies on their SDSS images are nearly compact 
at low redshifts and unresolved at high redshifts. Their 
typical angular sizes are less than 10". 



Izotov et al.l ([201 ll ) used all LCG spectra and Monte 
Carlo simulations to fit spectral energy distributions 
in the wavelength range AA3800 - 9200A. As for star 
formation history they assumed a single young burst 
with the age which is varied in the range < 10 Myr, 
and a continuous star formation with a constant star 
formation rate, which started at the lookback time t\ 
and finished at the lookback time ti < t±. Parame- 
ters t\ and t2 are varied in the range 10 Myr -13 Gyr. 
The contribution of gaseous continuum in LCGs is very 
high, therefore it was fitted first using equivalent width 
EW(H/3) of the H/3 emission line and subtracted from 
the observed spectrum prior fitting of the stellar con- 
tinuum. The masses of the young and old stellar pop- 
ulations, the age of the young burst and parameters 
for the old stellar population t\ and ti were parame- 
ters of fittin g . Mo re details of fitting can be found in 
Izotov et al.l (120111). We use t he results of modelling 



obtained bv llzotov et al.l (|2011l) . 

We use the GALEX Medium Imaging Survey (MIS) 
and All-sky Imaging Survey (AIS) data (see http: 
//galex. stsci.edu/GR4) to estimate the galaxy UV SFR 
for the LCGs sample. GALEX is a NASA Small Ex- 
plorer mission performed the all sky ultraviolet survey 
in two bands: far-UV (FU V, A e g= 1528A), and near- 
UV (NUV, A off = 2271 A) (|Morrissev et al.ll2005h . MIS 
and AIS data contain information on fluxes of ~10 7 
galaxies. The prime goal of GALEX is to study star 
formation in galaxies and its evolution with time. The 
major science objectives and c haracteristics of GA LEX, 
and of surveys are desc ribed bv lMartin et al.l (|2005f ) and 



Morrissev et all (|2005j). 

We matched the GAL EX data and the sample of 
LCGs (jlzotov et al.ll201lh and extracted FUV and NUV 
fluxes from the GALEX MIS and AIS database. These 
data combined with the NASA/IPAC Extragalactic 
Database (NED) data provide the determination of the 
galaxy UV luminosities. 

We excluded the galaxies with the UV flux errors 
exceeding 50% and the data for the multiple UV sources 
within the aperture of ~30". 



3 Correction for extinction 

Because the radiation of galaxies is reduced by dust ex- 
tinction, we applied r eddening corre c tions to Ha and 
UV band fluxes using ICardelli et aD (Il989h reddening 
law. Adopting the _R (U)-dependent extinction law from 
Cardelli et ail dl989t ) with R(V) = A(V)/E{B - V) 
= 3.1, we obtain A(Ha) = 2.54 x E[B - V) in Ha, 
A(FUV) = 8.l5xE(B - V) in the FUV band and 
4(NUV) = 9.17xE{B - V) in the NUV band. The 
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extinction coefficient C(H/3), reddening E(B — U)sdss, 
and the equivalent width of underl ying stellar hydrogen 
absorption lines were obtained bv llzotov et al.l ( 2011 ) 
from the hydrogen Balmer decrement in the redshift- 
corrected spectra. All hydrogen line fluxes were cor- 
rected for both the reddening and underlying stellar 
absorption. For comparison, we also use the redden- 
ing E(B — U)ned from the NED which was obtained 
from t he Milky Way reddening maps bv ISchlegel et al ' 



dl998h . The E(B -V) S nss and E(B-V) NE d differ be- 
cause the former quantity is the total reddening along 
the line of sight which includes extinction from both 
the Milky Way and the studied galaxy, while the latter 
q uantity is the r edden ing in the Milky Way only. 

llzotov et al.l (|201lh derived E{B - U) SD ss = for 
65 out of 803 galaxies. For some other galaxies they 
obtained E{B - U) S dss < E(B - U) N ed- To correct 
galaxy fluxes for extinction we use the E(B — U)sdss 
for all LCGs. Alternatively, we analyse all the UV data 
adopting E(B - V) SDS s if E(B - V) SDSS > E(B - 
U)ned and E(B - V)ned if E(B - U) S dss < E(B - 
U)ned- The difference between these two approaches is 
small because of the low extinction in LCGs and does 
not influence appreciably our results and conclusions 
(see Sect. g]). 

We derive mean reddenings E(B — U)sdss of 0.133 
and 0.134 for LCGs detected in the FUV and NUV 
ranges, respectively. We also find that the mean red- 
dening difference E(B — U)sdss — E(B — U)ned is 
0.106. This difference is the rough average internal ga- 
laxy reddening which can be used for galaxies without 
spectroscopic estimates of reddening. In principle, the 
correction for the Milky Way and intrinsic reddening 
should be done separately. This is because the Milky 
Way correction has to be applied to the fluxes at the 
observed wavelengths, while the correction for the in- 
trinsic reddening should be applied to the fluxes at the 
rcdshift-corrected wavelengths. However, ignoring of 
the separate correction for the Milky Way and intrinsic 
reddening would introduce very small additional uncer- 
tainties in the extinction-corrected fluxes. All LCGs are 
located at high galactic latitudes where the Milky Way 
extinction is very low, with the mean E(B — V)ned of ~ 
0.03. Most of LCGs are also low-redshift galaxies with z 

< 0.2. Only few galaxies are at redshifts z > 0.3. There- 
fore, the difference between the extinction correction of 
the flux applied with the observed wavelength and the 
redshift-corrected wavelength for the galaxy with the 
redshift z — 0.3 and adopting E(B — V)ned=0.03 is 

< 10% in the FUV and NUV bands, and <3% for the 
H/3 and Ha emission lines. This difference for galaxies 
with lower redshifts is lower, e.g. it is only ~ 5% in 
the FUV and NUV bands and - 2% for H/3 and Ha 



for the galaxy with z = 0.2. Furthermor e, Milky Way 
extinction maps bv ISchlegel et al.1 (|l998f ) are obtained 
for large apertures of 6', which are much larger than 
the angular sizes of LCGs. Th e refore , the extinction 
derived from the ISchlegel et al.1 (|l998l ) maps may not 
correspond to the real extinction in the direction on 
the galaxy if small-scale spatial extinction variations 
are present. These are reasons why we do not separate 
correction for the Milky Way and intrinsic reddening 
and use in a subsequent a nalysis E{B — U)sdss de- 
rived bv llzotov et alJ (|201ll ) from the hydrogen Balmer 



decrement. 

Accordingly to llzotov et al.l (|201lh we split our sam- 
ple into two subsamples of 276 "regular" galaxies with 
the round shape and 519 "irregular" galaxies with some 
sign of disturbed morphology suggesting the presence of 
two or more star-forming regions and their interaction. 
For these two subsamples we obtained E(B — U)sdss 
- E(B - V)ned of 0.081 and 0.120, respectively. Using 
the Student criterion we derived t = 4.98. This value 
suggests that the probability of the statistically signif- 
icant difference of E(B ~ U)sdss - E(B — U)ned for 
two subsamples is greater than 99.9%. Therefore, ex- 
tinction in galaxies with the non-round morphology is 
higher. On the other hand, we do not find tight cor- 
relation between E{B — U)sdss — E(B — U)ned and 
heavy element abundances. Apparently, the extinction 
is determined not only by the dust mass which is ex- 
pected to be higher in galaxies with higher metallicity, 
but also by the spatial distribution of dust. This distri- 
bution seems to be different in "regular" and "irregular" 
galaxies. 

We correct the galaxy fluxes for extinction accord- 
ing to I(X) = F(X) x 2.512 A( - X \ where F(X) and I(X) 
are the observed and the corrected fluxes, respectively. 
The extinction-corrected GALEX FUV and NUV fluxes 
from LCGs are nearly three times higher than the ob- 
served fluxes. In addition, Ha fluxes were corrected for 
an aperture comparing the total galaxy apparent mag- 
nitude m and the magnitude m(3") inside the SDSS 
spectroscopic aperture of 3" in a certain SDSS band 
depending on the galaxy redshift. We compare SDSS 
magnitudes m = r and m(3") = r(3") for galaxies with 
redshifts < 0.04, i and i(3") for galaxies with redshifts 
in the range 0.04 - 0.26, and z and z(3") for galaxies 
with redshifts > 0.26. Then, the aperture Ha flux cor- 
rection is A—2.512 m ^ m ^ > , where m — r,i,z depending 
on the galaxy redshift. 



4 Relations between galaxy luminosities and 
other global characteristics 

For each galaxy we calculated its Ha, FUV and NUV 
luminosities. We use the regression analysis to study 



Ha and UV luminosities and SFRs in luminous compact galaxies 



5 



a dependence of the LCG luminosities on other their 
characteristics. To provide a simple way for compar- 
ing L(Ha), L(FUV) and L(NUV) we use some other 
parameters being proportional to them. Namely, we 
use the calibration for SFRs averaged over the reason- 
able times c ale fo r different SFR tracers and defined by 
Kennicutt! jl998l ) as 



100 - 



SFR = txi. 



(1) 



The conversion factors k between the SFR and the 
L(Ra), L(FUV) and L(NUV) in Eq. [Dare derived us- 
ing the evolutionary synthesis models. The coefficient 
k depends on the time scale of star formation, initial 
mass function (IMF) and galaxy metallicity. Adopting 
the solar metallicity, the IMF with the powe r-law in- 
dex 2 .35 and mass limits of 0.1 and 100 Mq (jSalpeter 
Il955l) . SFR in M Q yr" 1 , L(Ha) in e r g s~\ L(FUV) and 
L(NUV) in erg s" 1 Hz- %ennicuttl (jl998l ) obtained the 
coefficient k of 7.9 x 10~ 42 for the Ha luminosity and 
1.4 x 10~ 28 for the FUV and NUV luminosities. 

Ther e are some other modific ations of Eq. [TJ In par- 
ticular, Kennicutt et al. ( 20091) proposed a composite 
SFR calibration based on the luminosities of both the 
Ha emi ssion line and t he UV continuum, adopting the 
IMF bv iKroupal f|200lh and obtained k = 5.5 x 10~ 42 . 
For clarity , we will use the values of k derived by 
Kennicutt ( 19981) for the solar metallicity and Salpeter 
IMF. Detailed review and analysis of SFR calibrations 
based on the Ha A6563A and [O II] A3727A emission 
lines, far infrared and ultraviolet con tinua are given by 
KennicuttJ jl998l ) and ICalzettJ (|2012l) . 

While the calibration defined by Eq. [TJ holds for 
continuous or quasi-continuous star formation which is 
common in the big galaxies with frequent starbursts, 
the situation is more complicated in the dwarf galaxies 
like LCGs with the strong and rare bursts of star forma- 
tion. In these systems with the instantaneous bursts the 
time interval of the star formation is not well defined, 
and the observed Ha luminosity strongly decreases on 
a time scale of a few Myr. Similar conclusions in lesser 
extent can be drawn for the SFRs derived from the 
FUV-band and the NUV-band luminosities. However, 
again for the clarity, we adopt that Eq. Q]can be applied 
for the determination of SFR in LCGs. 

A set of galaxy parameters includes the primary pa- 
rameters obtained directly from its SDSS spectrum, 
such as the redshift, the Ha flux and the chemical el- 
ement abundances. For each galaxy we also use sec- 
ondary parameters obtained by m odelling galaxy spec - 
tral energy distribution, following llzotov et al. I (|2oiil 
namely the masses of young M(young) and old M (old) 
stellar populations, the total mass of a stellar popula- 
tion M*, the age of a starburst i(young), and the lower 
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Fig. 1 Mass of the young stellar population m vs. the 
starburst age T. Dots and open circles correspond to sub- 
samples of "regular" and "irregular" galaxies, respectively 



(£2) and upper (ti) limits for the age of old stars. All 
these parameters are distance-independent. Neverthe- 
less, some nearby galaxies from our sample have angular 
diameters larger than the 3" SDSS aperture. Neglecting 
an aperture correction for these galaxies would result in 
the underestimation of stellar population masses. To 
estimate the proper masses we used the aperture cor- 
rection similar to the one used for luminosities. The 
luminosities of galaxies with angular diameters greater 
than 3" without this correction deviate substantion- 
ally from regression relations used for approximation of 
their SFRs as a function of their luminosities. These de- 
viations vanish after taking into account aperture cor- 
rections. 

We consider all above-mentioned parameters. How- 
ever, as one can see below, only two secondary param- 
eters, namely Al (young) and t(young) have the statis- 
tically significant impact on the luminosity. Hereafter 
we will use the brief notations m = M (young) and T = 
t (young) for the mass of the young stellar population 
and its age, respectively. 

Consider distributions of some secondary parame- 
ters. The values of the starburst age T are distributed 
rather uniformly over the interval T < 5 Myr. Galax- 
ies with T > 5.6 Myr are not included in our sam- 
ple due to the criterion EW(H/3) > 50A. This selection 
is also resulted in a relative decrease of the number 
of galaxies with 5 < T < 5.6 Myr. The distribution 
of m depends on T. This is illustrated in Figure [TJ 
One can see that the mean value of m for T > 3.2 
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Myr increases with the increase of T with the best 
fit log(m/lO 7 Af ) = -0.57 ± 0.22 + (0.33 ± 0.05) x 
T(Myr) for the subsample of "regular" galaxies and 
log(m/10 7 M Q ) = -0.75±0.17+(0.40±0.04)xT(Myr) 
for the subsample of "irregular" galaxies. 

Note that both values to and T were calculated as- 
suming a single star-forming region in the galaxy, while 
several regions of star formation with different m and 
T sometimes are observed in the galaxies. In this case 
we cannot rely on single values of to and especially T, 
the latter value would tend to be larger. To prove that 
we consider a case of the two star-forming regions with 
equal stellar masses, one is very young and another is 
older. The young star-forming region would dominate 
in the Ha luminosity because of the strong dependence 
of the flux of ionising radiation on a starburst age. On 
the other hand, the intensity of the optical continuum 
is less sensitive to the starburst age, therefore both 
star-forming regions equally contribute to the optical 
continuum, resulting in lower EW(Ha) as compared 
to the case with a single young burst. Consequently, 
this would result in a larger T because it is determined 
mainly by EW(Ha). It is natural to expect that galax- 
ies with larger masses of young stellar population with 
higher probability consist of several regions of star for- 
mation at different evolutionary stages. Therefore they 
would tend to have larger T as compared to the galax- 
ies with smaller masses of the young stellar population. 
That is why these galaxies concentrate in the upper 
right corner in Figure [T] 

We will show later that the ratio L/m decreases ex- 
ponentially with increasing T if T > 3.2 Myr. As a 
result, the lower right corner in Figure [1] is empty be- 
cause the sample is flux-limited L(H/3) > 3xl0 40 erg 
s _1 . Galaxies with low m and T > 4 Myr have lu- 
minosities below the threshold and do not enter the 
sample. 

Our goal is to find simple but statistically significant 
dependences of the galaxy luminosity on other primary 
and/or secondary galaxy parameters. First, we search 
for a set of parameters to which the galaxy luminosity is 
the most sensitive. Later, we will find the best formulae 
to describe these dependences and analyse them. 

At first we do not take into account galaxy metal- 
licities and analyse linear dependences of their lumi- 
nosities on other parameters and their combinations. 
Regressors were chosen accordingly to the statistical 
Fisher test (|Fisherlll954T) . We reject the regressors with 
statistical significance below the threshold of 99.95% 
and consider the regressions which are good for all six 
subsamples by selecting 3 wavelengths (Ha, FUV and 
NUV) and 2 morphologies - "regular" and "irregular" . 
These regressions have the form 



SFR = Ci + C 2 m + C 3 mT 2 + C 4 m 2 . 



(2) 



The values, the errors and the statistical significances of 
the coefficients Cj obtained by the least squares method 
(LSM) are shown in Table Q] for to expressed in solar 
masses and T in yr. The root mean square (RMS) 
standard deviation a and the number N of the galaxies 
in each subsample are also shown in the Table. Sig- 
nificances are characterised by the value F obtained 
by Fisher's test. The critical F values, corresponding 
to the statistical significances of 90, 95, 97.5, 99, 99.5, 
99.9 and 99.95% are equal to 2.71, 3.84, 5.02, 6.64, 7.88, 
10.83 and 12.10, respectively. One can see from Table[T] 
that the significances of regressors Ci and C3 in Eq. 
are higher than 99.95%. For the regressor C4 we choose 
the threshold value F > 10. Only for one subsample 
the threshold attains a higher value. For all other sub- 
samples we assume C4 = and indicate in parentheses 
the value of F for the case C4 ^ 0. 

Consider an implication of Eq. [2] One would expect 
that SFRs vanish at the low-mass limit to = 0. How- 
ever, Eq. [2] (case (a) in Table [1]) implies that SFR is 
equal to non-zero C\ at this limit with a large statistical 
significance for the three subsamples out of six. We as- 
sume that this is due to the uncertainties in the estima- 
tion of m. We suggest that this statistical effect resem- 
bles the well-known Malmquist bias. A similar effect 
was studied in con nection with the large-scale co llec- 
tive galaxy motion (jParnovskv fc Parnowskill2005l ). To 
verify this hypothesis we performed some Monte Carlo 
simulations. For these simulations we need many gener- 
ated mock catalogues, preferably with the distribution 
of the parameters similar to that in real subsamples. 

First, we adopt the values of to and T from the cor- 
responding real subsamples. Then, we set G\ = and 
calculate the coefficients C2 , C3 , and C4 (if the last one 
is statistically significant) for the regression relation Eq. 
[2] using the LSM. These coefficients are shown in Table 
[U case (b). After that, we derive the SFR values for 
the case (b) from Eq. [2] with C\ =0. As a result we 
obtain the initial "unbiased" set of to, T and SFR val- 
ues. Note that the dependence of SFR on m and T in 
this sample is functional, not statistical. The real SFR 
data are taken into account only indirectly via the set 
of coefficients C2, C3 and C4. Naturally, if we apply 
the LSM with the regression defined by Eq. [3] to this 
data we obtain coefficients in the corresponding row of 
Table [U case (b), but with a = 0. 

Next, using the Monte Carlo technique we add a 
noise, i.e. random errors to the unbiased values of SFR 
or to. We find that random normal errors in SFR val- 
ues result in the nonshifted distributions of Cj values 
obtained by the LSM. In this case a random value for 
C\ has the low statistical significance. A completely 
different situation arises when random errors in m are 
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Table 1 Coefficients C, in Eq. [5] and their significance values F for different galaxy subsamples 



Subsample 


N 


a 


Ci(F) 


C* 2 x 10 8 (F) 


C 3 x 10 22 (F) 


C 4 x lO 18 ^) 


a) the case with C\ ^ 


1. 


Ha, "regular" 


276 


3.4 


1.15 ±0.35(10.9) 


18.3 ±0.8(528) 


-43.1 ± 2.6(268) 


-33.8 ±8.0(17.8) 


2. 


Ha, "irregular" 


519 


4.0 


1.41 ± 0.24(35.3) 


17.0 ±0.4(1602) 


-45.5 ± 1.6(766) 


-(0.7) 


3. 


FUV, "regular" 


213 


4.8 


0.19 ±0.45(0.2) 


12.4 ± 1.0(149) 


-22.1 ±4.0(30) 


-(3.9) 


4. 


FUV, "irregular" 


418 


4.4 


1.29 ±0.31(17.5) 


8.0 ±0.6(159) 


-14.2 ±2.4(34) 


-(0.4) 


5. 


NUV, "regular" 


233 


5.6 


0.42 ±0.50(0.7) 


14.5 ± 1.1(167) 


-20.7 ±4.4(22) 


-(2.3) 


6. 


NUV, "irregular" 


435 


7.8 


0.27 ±0.53(0.3) 


14.8 ± 1.1(180) 


-23.5 ±4.2(31) 


-(8.7) 


b) the case with G\ = 


1. 


Ha, "regular" 


276 


3.5 





20.2 ± 0.6(1275) 


-46.3 ± 2.5(343) 


-49.2 ±6.6(55.2) 


2. 


Ha, "irregular" 


519 


4.1 





18.4 ± 0.4(2592) 


-49.8 ± 1.5(1062) 


-(2.0) 


3. 


FUV, "regular" 


213 


4.8 





12.6 ±0.8(227) 


-22.8 ±3.6(40) 


-(1.8) 


4. 


FUV, "irregular" 


418 


4.5 





9.6 ± 0.5(352) 


-19.3 ±2.1(82) 


-(0.4) 


■5. 


NUV, "regular" 


233 


5.6 





15.1 ±0.9(264) 


-22.3 ±4.0(32) 


-(0.5) 


6. 


NUV, "irregular" 


435 


7.8 





15.2 ±0.9(303) 


-24.6 ±3.6(47) 


-(7.1) 



considered. In this case the distributions of Ci values 
are shifted relative to the "unbiased" ones and we ob- 
tain a non-zero value for C\ , sometimes with the large 
false statistical significance. Other coefficients tend to 
attain values nearer to zero if the noise increases. 

Could the effect of uncertainties in the to determi- 
nation explain a non-zero value of C\ obtained from 
the real data? To prove this suggestion we compare 
the values of Ci in Table [1] (case (a)) with the ones 
obtained using Monte Carlo simulations. For that, we 
use the "unbiased" SFR and T values from the sample 
calculated with the coefficients from case (b) in Table [1] 
and add a random noise to to. We obtain the "biased" 
value rribias- The distribution of muas is expected to 
be log-normal, therefore we added the noise using the 
equation 

muas = m(l + s x f), (3) 

where £ is a normally distributed random value with 
zero mean and unity dispersion and s characterises an 
amplitude of the noise. Then we calculate Ci values ap- 
plying LSM to the set of "unbiased" SFR and T values 
and the values of rribias- After repeating this procedure 
1000 times we obtain the mean values and the distribu- 
tion of Ci as well as the mean value of a. The results 
for different s are shown in Table [5] for the Ha subsam- 
ple No. 1 in Table [TJ One can see that the calculated 
values of C\ in Table Q] (case (a)) can be explained by 
the Monte Carlo model with s = 0.34 (Table [2). In 
addition to the shift of C\ the noise in to also affects 
the values of other coefficients Cj. They are shifted 
closer to the values from Table Q] (case (a)). The value 
of modelled a is smaller than that obtained from real 
data due to the contribution of SFR errors in real data, 
in addition to errors in to. 



We used similar Monte Carlo simulations for other 
two subsamples with statistically significant non-zero 
C\ values (subsamples No. 2 and 4 in Table [1} and 
obtained similar results attained at s = 0.32 and s — 
0.38, respectively. 

Therefore, we adopt that the true preliminary re- 
gression for SFR is Eq. [2] with C\ = and non-zero 
values of Ci, C3 and sometimes C4 (see case (b) in Ta- 
ble [T]) . The uncertainties in the estimation of to lead 
to the appearance of the first regressor in Eq. [2] with a 
false statistical significance. Hereafter we will use only 
the models with C\ = 0. Then we can rewrite Eq. [5] in 
the form 

SFR/to = C 2 + C3T 2 + C 4 m, (4) 

introducing a new value SFR/to which gives us a pos- 
sibility to consider the dependence on T regardless of 
the dependence on to. The first term in Eq. |4] is the 
main one and its meaning is that SFR of a galaxy is 
approximately proportional to the mass of the young 
stellar population. This is quite obvious because ra- 
diation in the UV continuum and Ha emission line is 
emitted mainly by young O-stars. We will return later 
to the consideration of the possible nonlinear depen- 
dence of SFR on to. 

In Figure [2] we show the dependence of L(Ha)/TO 
oc SFR/m on the starburst age T. It is seen from the 
Figure that the ratio L(Ha)/m is practically constant 
for T < To =3. 2 Myr and decreases practically exponen- 
tially for larger T. We note that SFR(Ha), to, and T 
are not directly correlated because they are based on 
the different features in the spectra: SFR(Ha) is not 
modelled and is derived from the Ha flux, to and T are 
modelled but they depend on different features. The 
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Table 2 Coefficients d in Eq. [5] for the subsample No. 1 (Table [TJ obtained from Monte Carlo simulations adopting 
"unbiased" values of SFR and T and different amplitudes s of the noise for m according to Eq. [3j 



s 


Ci 


C 2 x 10 8 


C 3 x 10 22 


d x 10 18 


17 


0.20 


0.44 


19.1 


-43.1 


-45.3 


1.83 


0.30 


0.98 


17.5 


-39.0 


-40.9 


2.72 


0.32 


1.11 


17.2 


-38.1 


-39.8 


2.87 


0.34 


1.24 


16.9 


-37.2 


-38.7 


3.03 


0.35 


1.30 


16.6 


-36.5 


-38.2 


3.10 


0.40 


1.61 


15.8 


-34.3 


-35.4 


3.43 



E 

S 
X 

DC 

w 




Fig. 2 Ratio L(Ha)/m of the galaxy luminosity in the Ha 
emission line to the mass m of the young stellar popula- 
tion vs. the age of the starburst T. Dots and open circles 
correspond to subsamples of galaxies with "regular" and 
"irregular" shape. The solid and dashed lines correspond 
to the best fits for the subsamples 1 and 2 with the mean 
galaxy metallicity using the regression Eq. [9] 



mass of young stellar population to is determined by 
the luminosity and the shape of the stellar continuum 
in the visible range, while T depends mainly on the 
equivalent width EW(H/3) of the H/3 emission line. In 
particular, the equivalent width EW(H/3) is lower by 
a factor of many times for the burst with the age of 6 
Myr as compared to that for the burst with the age of 3 
Myr, while the flux of th e stellar continuum ne ar H/3 is 
decreased by only ~ 5% (jLeitherer et al.lll999l ). There- 
fore, SFR(Ha), m, and T can be considered as indepen- 
dent parameters. The most plausible explanation of the 
L(Ha)/m - T relation is that most massive stars with 
masses of ~ 100 Af Q , producing most of the ionising ra- 
diation, disappear after the starburst age T ~ 3.2 Myr, 
corresponding to their lifetimes. It also implies that 
stars in LCGs are formed during very short time peri- 



ods, otherwise, in the case of extended bursts, To would 
be greater than ~ 3.2 Myr. The best fit for T > 3.2 
Myr is log(SFR/m) = -5.61±0.07+(-0.316±0.016) x 
T(Myr) for the "regular" subsample and log(SFR/m) = 
-5.62 ± 0.05 + (-0.326 ± 0.011) x T(Myr) for the "ir- 
regular" subsample. 

Is this effect statistically significant? The values 
of the Fisher coefficients for linear terms in regression 
relations for both subsamples exceed 398. For addi- 
tional proof we performed the Student test to com- 
pare the mean values of the SFR/to for T < 3.5 Myr 
and T > 4 Myr. In all cases the Student test gives 
t above 3.31, corresponding to the statistical signifi- 
cance of 99.95%. The results are shown in Table [31 
where a is the standard deviation. One can see that 
mean values of SFR/to for T < 3.5 Myr are essen- 
tially larger than those for T > 4 Myr. This effect 
is more pronounced for the Ha radiation as compared 
to the FUV and NUV radiation. This is because 1) 
more massive short-lived stars contribute to the ionis- 
ing radiation and 2) the luminosity of ionising radiation 
is much stronger increased with the mass of a star as 
compared to the UV-radiation. We note, however, that 
the ratio SFR/m decreases more slowly in comparison 
with the Schaerer fc Vacca ( 1998t) population synthesis 
models for young stellar populations given the appro- 
priate heavy element abundance. 

In accordance with the above discussion we introduce 
the function f(T) 



/CO 



1 



exp(-p x (T - T )) 



if T<T ; 
if T>T , 



(5) 



where To = 3.2 Myr and p = 0.75 Myr^ 1 . These values 
are the preliminary ones obtained from Figure [5] They 
will be improved later. Using f(T) and transforming 
Eq. 0] we obtain the regression relation 



SFR = C 5 mf(T) + dm 2 



(6) 



and apply LSM to calculate the values and the errors 
of the coefficients. 
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Table 3 Comparison of the average values of SFR/m for different starburst ages T according to the Student's t-test 



Sample 




T < 3.5 Myr 






T > 4 Myr 




t 


N 


SFR/m x 10 s 


a x 10 s 


N 


SFR/m x 10 s 


a x 10 s 


1. Ha, "regular" 


137 


20.9 


6.9 


100 


8.8 


3.3 


17.8 


2. Ha, "irregular" 


249 


20.9 


8.1 


215 


7.9 


3.0 


23.5 


3. FUV, "regular" 


104 


12.4 


8.7 


79 


7.5 


5.9 


4.6 


4. FUV, "irregular" 


198 


13.8 


11.8 


177 


5.4 


2.5 


9.7 


5. NUV, "regular" 


114 


15.6 


11.1 


88 


9.7 


7.5 


4.5 


6. NUV, "irregular" 


207 


19.5 


18.6 


181 


8.4 


4.9 


8.2 



The improved dependence of SFR/m on T (dashed 
and solid lines in Fig. [2] correspond to best fits with the 
optimal values of p) and the existence of the correlation 
between T and m (Fig. [TJ lead us to the following re- 
sults. For the Ha "regular" subsample the second term 
becomes statistically insignificant and for Ha "irregu- 
lar" subsample the second term becomes positive and 
statistically significant. The values of RMS standard 
deviations become smaller than those from Eq. [4] be- 
cause of the better approximation of SFR/m on T. 

The dependences of L(FUV)/m and L(NUV)/m on 
T are more gently sloping and have a larger scattering of 
L(FUV)/m and L(NUV)/m as compared to L(Ha)/m. 
The dependences L(FUV)/m on T and L(NUV)/m on 
T are shown in Figure |3l The shape and application 
of f(T) to describe these dependences will be discussed 
below. First of all we will use the function f(T) for the 
Ha radiation. 

We now consider the dependence of L on the galaxy 
metallicity. As a measure of metallicity we choose the 
oxygen abundan ce [O] = 12 + log(0 /H). [O] was accu- 
rately derived bv llzotov et al. I ((20li for every galaxy in 
the LCG sample. We adopt these values for the regres- 
sion determination. To analyse the trend we restrict 
ourselves to the linear dependence on [O] only. Oxygen 
abundances of the galaxies from the LCG sample vary 
in the range from 7.52 to 8.47 with the median value of 
8.06. We note that the use of nonlinear dependences on 
[O], for example, 10^ , does not make any additional 
progress in reducing of the regression RMS standard 
deviations. 

First, we discuss potentially misleading methods of 
studying the trends in the dependence of L on [O] . The 
simplest approach is to calculate the mean L for the 
subsamples with low and high oxygen abundances and 
to treat the obtained difference as an dependence on 
the oxygen abundance. Adopting the median value [O] 
= 8.06 as a value dividing galaxies with low and high 
oxygen abundances, we find that the subsamples with 
high [O] have the mean L (or SFR) values ~ 1.5 times 
larger than for the subsamples with low [O]. However, 
the mean masses of young stellar population in higher- 
metallicity subsamples are larger by a factor of ~ 1.5 



in comparison with subsamples with lower metallicity. 
Therefore, the differences in SFRs are mostly due to the 
differences in masses of the young stellar population 
and are not directly related to the differences in the 
metallicity. 

The best way to investigate the direct impact of ga- 
laxy metallicity on its luminosity for the galaxy with 
the same values of m and T is to include the metallic- 
ity directly in the regression relation. Therefore, we fit 
the SFR values for Ha emission line using the set of 
regressors 



SFR = C 5 mf(T) + dm 2 + C 6 m([0] - ([O])) 



(7) 



with T = 3.2 Myr and p = 0.75 Myr" 1 . Here ([O]) 
is the mean oxygen abundance of the sample galax- 
ies. We subtract this value from the galaxy's oxygen 
abundance [O] to make the last term in Eq. [7] prac- 
tically orthogonal to the first one and in this way to 
keep the results of the Fisher test for the first two re- 
gressors. Using the LSM we calculate the coefficients, 
their errors and statistical significances of the regres- 
sion defined by Eq. [7J For the Ha "regular" subsample 
we obtain a = 3.42, C 5 = (2.13 ± 0.05) x 10~ 7 (F = 
1654), C 4 = (-3.6 ± 4.2) x lO^ 18 ^ = 0.7), C 6 = 
(-4.7 ± 1.0) x 10~ 8 (F = 20). For the Ha "irregu- 
lar" subsample the derived values are a = 3.8, C5 = 
(1.94 ± 0.03) x 10~ 7 (F = 4415), C 4 = (3.4 ± 1.6) x 
10~ 18 (F = 4.9), C* 6 = (-3.2 ±0.7) x lO" 8 ^ = 23). We 
note the drop of statistical significance of the nonlinear 
term with C4 below the threshold for both subsamples. 

Switching to the UV luminosities, we generalise Eq. 
[4] by adding the term with the dependence on metallic- 
ity: 

SFR = C 2 m + C 3 mT 2 + C 4 m 2 + C 6 m([0] - ([O])). (8) 

Derived coefficients, their errors and statistical signif- 
icances for FUV and NUV subsamples are shown in 
Table |4j All coefficients C§ are negative and have the 
statistical significanse more than 99.5%. Coefficients 
C4 are statistically insignificant for all FUV and NUV 
subsamples, similar to that for Ha subsamples. 
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Fig. 3 Ratio L(UV)/m of the luminosity in the UV range to the mass m of the young stellar population vs. the age of 
the starburst T for the a) FUV and b) NUV radiation. Dots and open circles correspond to subsamples of galaxies with 
"regular" and "irregular" shape. The solid and dashed lines show the best fits (Eq. [9} with f(T) in Eq. [5] and Eq. I10I 
respectively, and adopting [O] = ([O]) 



Summarising, we find that a nonlinear term m 2 in 
the regression relations (Eqs. [7] and [8]) is statistically 
insignificant for all six subsamples. Speaking more pre- 
cisely, we conclude that the statistical analysis gives us 
no reason to justify the existence of such a term. Would 
it be statistically significant it makes the ratio L/m be 
depending on the young stellar population masses. This 
could be treated as an impact of the environment or as a 
result of some kind of an interaction of several regions of 
star formation. However, it is difficult to analyse these 
effects statistically because of the correlation between 
m and T, which would lead to ambiguous conclusions. 

We discard nonlinear term with m 2 adopting C4 
I ) in Eq. [7] and obtain a regression relation in the form 

SFR = C 5 mf(T) + C 6 m{[0}-([0]}). (9) 

We apply this relation for the Ha radiation. For the 
subsamples No.l and 2 we use the function /(T) from 
Eq. with T = 3.2 Myr and p = 0.75 Myr" 1 . The val- 
ues, errors and the statistical significances of the coeffi- 
cients obtained by the LSM are shown in Table [5] (case 
(a)). Using these coefficients we calculate SFR resr (Ha) 
and L regr (Ra) = SFR regr (Ha)//c for every galaxy from 
these subsamples, where k = 7.9 x 10~ 42 . The compar- 
ison of the calculated values with the observed ones is 
plotted in Figure |3J It follows from the Figure that Eq. 
O provides a good approximation of the observational 
Ha data in the entire range of SFR(Ha) =0.8-77 M 
yr" 1 . 

For the FUV and NUV subsamples we use Eq. [5] 
with C4 = 0. We can rewrite it in the form of Eq. [SJ 



introducing 

f(T) = I - V T 2 , C b = C 2 , V = (10) 

From the two last rows of Table [4] (case (b)) we derive 
77 = (21.6±0.9) x 10~ 3 Myr^ 2 for the FUV band and 77 = 
(18.5±0.9) x 10~ 3 Myr~ 2 for the NUV band. Certainly, 
the relation Eq. [TUlcannot be used if the starburst age T 
is greater than 7 Myr because f(T) becomes negative. 
Note that the errors of 77 were estimated taking into 
account not only the errors of C2 and C3 in Table 0] 
but their covariation too by using all elements of the 
correlation matrix including the nondiagonal ones. 

Consider the final fine tuning of the parameters in 
the Eq. [9] Using LSM we obtain the optimal values 
of To and p in Eq. [5] For the Ha "regular" subsample 
we derive a — 3.3, T = 3.3 Myr, p = 0.69, while for 
the Ha "irregular" subsample the derived values are 
a = 3.6, T = 3.3 Myr, p = 0.68. 

Could Eq. [9]be used for the UV continuum radiation 
with /(T) in the form of Eq. [5l similar to the Ha 
radiation? Such attempt turns out to be successful. For 
the subsamples No. 7 and 8 we obtain the optimal values 
T = 2.9 Myr and T = 3.1 Myr, respectively. They are 
smaller than 3.2 Myr, probably, due to the larger data 
scatter in comparison with the Ha subsamples. With 
these optimal values of To and p, the RMS standard 
deviations for the UV bands are slightly decreased as 
compared to the case when /(T) is used in the form of 
Eq. |TD] The minimum of the sum of square residuals 
^(SFRj — SFR) 2 , corresponding to the best value of 
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Table 4 Values of the coefficients in Eq. [H] 



Subsample 


N 


a 


C 2 x 10 s (F) 


C 3 x 10 22 (T) 


c 4 


x 10 18 (F) 


C 6 x 10 8 (T) 


([O]) 








a) Rc£ 


;ressions with C4 


^0 








3.FUV,regular 


213 


4.7 


13.6±0.9(231) 


-31.9±4.1( 60) 


30.9 ± 10.5(8.6) 


-7.0 ±1.7(16.2) 


8.05 


4. FUV, irregular 


418 


4.4 


10. 2±0. 5(375) 


-20.7±2.7( 60) 


-3.4 ±3.5(1.0) 


-3.9 ± 1.0(16.7) 


8.13 


5. NUV, regular 


233 


5.5 


16.0±1.0(246) 


-30.2±4.7( 42) 


23.8 


± 12.1(3.9) 


-5.8 ±1.9(9.1) 


8.05 


6. NUV, irregular 


435 


7.3 


17.8±0.9(405) 


-39.1±4.3( 83) 


12.9 ±5.6(5.3) 


-12.0 ± 1.6(57.7) 


8.13 








b) Rej 


;ressions with C4 


= 








3.FUV,regular 


213 


4.7 


13.8±0.9(231) 


-26.7±3.8( 50) 







-4.9 ± 1.6(9.2) 


8.05 


4.FUV,irregular 


418 


4.4 


10.3±0.5(378) 


-22.3±2. 2(101) 







-3.8 ±0.9(16.1) 


8.13 


5.NUV,regular 


233 


5.5 


16.1±1. 0(248) 


-26.2±4.2( 38) 







-4.2 ± 1.7(5.7) 


8.05 


6 . NUV ,irregular 


435 


7.3 


17.3±0.9(403) 


-33.5±3.6( 88) 







-12.2 ± 1.6(59.9) 


8.13 


7.FUV,all 


631 


4.8 


11.7±0.5(575) 


-25.2±2.0(157) 







-4.8 ±0.9(31.7) 


8.10 


8.NUV,all 


668 


6.8 


17.1±0.7(647) 


-31.6±2.8(130) 







-9.7 ± 1.2(65.4) 


8.10 



Table 5 Values of the coefficients in Eq. [9] 



Subsample 


N 


er 


P, Myr- 1 


C 5 x 10 7 (T) 


C 6 x 10 8 (T) 


([0]) 






a) Re 


gressions with p — 0.75 Myr 1 






l.Ha, regular 


276 


3.4 


0.75 


2.10±0.04(2752) 


-4.8±1.0(21) 


8.05 


2. Ha, irregular 


519 


3.8 


0.75 


1.98±0.02(7156) 


-3.4±0.7(27) 


8.13 






b) Regressions 


with the optimal p 






l.Ha, regular 


276 


3.3 


0.66 


1.95±0.04(2913) 


-4.6±1.0(21) 


8.05 


2. Ha, irregular 


519 


3.6 


0.65 


1.82±0.02(7849) 


-3.2±0.6(26) 


8.13 


9.Ha, all 


795 


3.5 


0.65 


1.85±0.02(10851) 


-3.6±0.5(46) 


8.10 


7. FUV, all 


631 


4.7 


0.42 


1.11±0.02(2771) 


-5.1±0.8(39) 


8.10 


8. NUV, all 


668 


6.7 


0.32 


1.60±0.03(3925) 


-10.1±1.1(79) 


8.10 


c) Regressions with the optimal p and alternative correction for extinction 


7. FUV, all 


631 


4.7 


0.43 


1.12±0.02(2799) 


-5.2±0.8(41) 


8.10 


8. NUV, all 


668 


6.7 


0.33 


1.63±0.03(3981) 


-10.3±1.1(83) 


8.10 



To, is rather shallow. Therefore, we can adopt a single 
value T = 3.2 Myr for all subsamples. 

In Table [5] (case (b)) we show the final values of the 
parameters for the subsamples No. 1, 2, 7, 8 and 9 in 
the form according to Eqs. [9] and [5] 

One can see that the approximation Eq. [5] is much 
better than Eq. [TU] for the FUV band and is slightly 
better for the NUV band. Moreover, the use of Eq. [5] 
for f(T) in Eq. [9] is more preferable not only because 
of the decrease of the RMS standard deviation a, but 
also by the same dependence on T as that in the case of 
Ha radiation. Starting from To, the fading half-times 
of the Ha, FUV and NUV emission are 1.1 Myr, 1.6 
Myr and 2.1 Myr, respectively. 

It is noted in Sect. [3] that for the reddening correc- 
tion of galaxy fluxes we use E{B — U)sdss derived from 
the hydrogen Balmer decrement. Alternatively, we also 
consider reddening corrections, adopting E(B — U)sdss 
if£(.B-U)sDss > E(B-V)kev (for -90% of the sam- 
ple) andE(B-V) NE D X E(B- V) S vss < E(B-V) NEB 



(~ 10% of the sample). The coefficients in this case are 
shown in Table (case (c)). The comparison of case 
(b) and case (c) coefficients shows that differences are 
very small, indicating that both approaches can equally 
be used. 



5 Star formation rates 

We already noted in Section [4] that there are different 
indicators of star formation in a wide range of wave- 
lengths from UV to radio and different calibrations to 
quantify it. As it is mentioned in numerous papers, 
each SFR indicator possesses its own strengths and 
disadvantages. Recently, the hybrid SFR indices were 
proposed, which are based on the combination of the 
ultraviolet and infrared tracers, the Ha and the in- 
frared or radio continuum tracers, the [O 11] A3727A 
forbidden-line doublet and the infrared or radio contin- 
uum tracers. Studies of star formation rates for differ- 
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Fig. 4 Luminosities in the Ha line obtained from the re- 
gression relation Eq. [9] with p = 0.75 Myr - . vs. the mea- 
sured ones. Dots and open circles correspond to subsamples 
of galaxies with "regular" and "irregular" shape 



ent samples of galaxies with different level of star for 
mation activity and with different SF tracers were car 



ried out in many papers (see, e.g. 



Boauien et all 2009; 



Calzetti et al.H2010t iGilbank et al ll2010l; iHopkins et al. 
2002; Iglcsias- Paramo e t al. 2006; Izotov a fc Parnovskv 



2006; iKennicutt ll998tlKennicutt et al.|2009l:lKewlev et alM 



2009; iLee et al.ll2009t iLi et al.ll20ld Moustakas i et al 



2006; Overzier et al. 2009; Sargsvan fc Wccdman 2009; 



Schmitt et alj|2006h . 

In the present paper, the galaxy luminosities in the 
Ha emission line and in the UV non-ionising contin- 
uum are used to obtain SFRs (Eq. [T]). The Ha emis- 
sion in the star-forming galaxies is produced by the gas 
ionised by the most massive short-lived hot O-stars with 
masses > 17 M and traces the star formation over the 
period of a few Myr, corresponding to the lifetime of 
these stars. The non-ionising UV emission is produced 
by stars in a wider range of masses and therefore can 
in principle be used as a SFR tracer on a time scale 
of up to 100 Myr. However, in the case of strongly 
star-forming LCGs, similar trends in Figs. [2] and [3] 
imply that Ha, FUV and NUV emission in LCGs are 
produced by the same young stellar populations. This 
conclusion is supported by the fact that instantaneous 
burst with the age of 6 Myr emits ~ 4 times and ~ 3 
times less radiation in the FUV and NUV ranges, re- 
spectively, as compared to that in the burst with the 
age of 3 Myr (jLeitherer et al.lll999l ). Similar difference 
is seen in Fig. [3] Calculating SFRs from Eq. Q]we ac- 



1 10 

SFR(Ha) 

Fig. 5 Current luminosities and SFRs in the FUV range 
vs. luminosities and SFRs in the Ha emission line. Dots 
and open circles correspond to subsamples of "regular" and 
"irregular" galaxies, respectively. Solid line is the line of 
equal SFRs, dashed lines show the best fits for subsamples 



tually use the galaxy luminosities observed at a certain 
current moment. However, in star-bursting galaxies, 
the observed Ha and UV-luminosities depend on the 
burst age and may vary on a time scale of several Myr. 
'This effect is most pronounced for the Ha luminosity: 
it is constant over first ~ 3 Myr of a starburst and then 
quickly declines with time. To take into account the 
temporal luminosity evolution we introduce the initial 
value of the Ha luminosity Lo(Ha) = L(Ha)(T — 0) af- 
ter the onset of star formation and calculate the value 
of SFR (Ha) = k x L (Ha). 

Similarly, we also introduce the initial luminosi- 
ties Lo(FUV) and L (NUV) in the FUV and NUV 
ranges and the corresponding values SFR (FUV) and 
SFRq (NUV) according to Eq. [1] Hereafter we consider 
the temporal evolution in UV ranges in the form of Eq. 
[3 To distinguish the functions f(T) for the Ha line 
and the FUV and NUV ranges we will use the corre- 
sponding subscripts. These functions differ only in the 
terms of the coefficient p values, which are presented in 
the Table El case (b). As an illustrative example, we 
will demonstrate now the certain advantage of using 
SFRo(Ha) in comparison with SFR(Ha). 

We have an ample sample of the galaxies with known 
Ha and FUV luminosities. First, we derive SFRs from 
the observed luminosities. We show in Figure[5]the rela- 
tion between SFR(FUV) and SFR(Ha) (or equivalently 
L(FUV) and L(Ha)). It is seen that both SFR(Ha) 
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Fig. 6 Initial luminosities and SFRo in the FUV range 
at zero starburst age vs. initial luminosities and SFRo in 
the Ha emission line. Dots and open circles correspond to 
subsamples of "regular" and "irregular" galaxies, respec- 
tively. Upper solid line is the line of equal initial SFRs, 
dashed lines show the best fits for subsamples, lower solid 
line SFR (FUV)= 0.52xSFR (Ha) shows their mean ratio 



and SFR(FUV) are of the same orders. The values of 
the LSM slopes in the dependences log L(FUV) on log 
-L(Ha) for subsamples No. 3 and 4 are 1.07 ± 0.05 and 
1.04 ± 0.04 (see dashed lines in Figure [5]). However, 
the values of the slopes in the inverse dependences log 
L(Ka) on log L(FUV) are much smaller than unity, 
0.66 ±0.03 and 0.72 ±0.03 for subsamples No. 3 and 4, 
respectively. These values would correspond to the lines 
with the slopes 1/0.66 = 1.51 and 1/0.72 = 1.40 in Fig- 
ure [5l The differences in the slopes of the direct and 
inverse dependences are likely due to the Malmquist 
bias caused by our selection of only galaxies with high 
L(H/3) > 3 x 10 40 erg s _1 . 

We reduce Ls and SFRs to a zero age of a starburst 
according to 



L = Lj f{T), SFRo = SFR//(T) 



(11) 



with /ff Q (T) for Ha emission line and /fuv(T') for 
FUV radiation. Adopting the values of p from Table 
[5j case b) , we obtain the distribution of galaxies shown 
in Figure [H The LSM gives the slopes 1.04 ± 0.04 and 
0.95±0.03for "regular" and "irregular" galaxies respec- 
tively, implying that the initial galaxy FUV luminosity 
Lo(FUV) is proportional to the initial Ha luminosity 
Lo(Ha). Slopes of the inverse linear dependences log 
SFR (Ha) on log SFR (FUV) for the same subsamples 



are 0.76 ±0.03 and 0.85 ±0.03. These values correspond 
to the lines with the slopes 1/0.76 = 1.31 and 1/0.85 = 
1.18 in Figure [HI Thus, though Malmquist bias is also 
present for the data reduced to the zero starburst age, 
its effect is much smaller because of smaller differences 
between the slopes of the direct log SFRo (FUV) - log 
SFR (Ha) and inverse log SFR (Ha) - log SFR (FUV) 
dependences. We also note that the data point scatter 
in Figure [6] is slightly smaller than that in Figure [5j 

However, there is a downward shift of the data points 
relative to the line of equal SFRs in Figure [6] indicating 
that SFRs obtained from different indicators are pro- 
po rtional, but not equ al . 

Ilzotov et al.l ( 2011 ) derived a single value of redden- 
ing E{B — V) for both the gaseous and stellar emis- 
sion assuming uniform distribution of dust. However, 
dust in galax i es is distri buted non-unifo r mly. E.g. 
Calzetti et all |l994h and ICharlot fc Fall (|2000l ) sug- 
gested that young massive stars responsible for the Ha 
emission are located in more dusty regions as com- 
pared to the stars which produce non-ionising UV radi- 
ation, including FUV and NUV ranges. In particular, 
Calzetti et al.l (|1994T) assumed that non-ionising UV ra- 
diation is produced by older stars which were formed in 
regions different from those where most massive young 
stars are present. 

Could the non-uniform distribution of dust explain 
the downward shift in Fig. [5]? Apparently, not. As- 
suming that E{B — V) for FUV and NUV ranges is 
smaller we obtain lower FUV and NUV luminosities. 
Therefore, the downward shift would be larger. Fur- 
thermore, as it was already noted above, similar trends 
in Figs. [2] and [3] imply that Ha, FUV and NUV emis- 
sion in LCGs are produced by the sam e young stellar 



popul ations, contrary to assumption by ICalzetti et al 
(|l994h . 



In order to equalise the values of SFR (Ha) and 
SFR (FUV) the coefficients k in Eq. Q] should be 
changed from their nominal values. We find that the 
mean ratios of SFR (Ha) /SFR (FUV) and SFR (Ha) 
/SFRo(NUV) are equal to 1.9 and 1.5, respectively. 
More precisely, we find that ioMsfRo(fuv)/sfr o( h q ))) 
= 0.53, io( 1o s( sfr o( nuv )/sfr (h q ))) = q gg 

Using these average ratios we obtain statistical rela- 
tions 



Lo(FUV) = 3.0 x 10~ 14 x L (Ha), 
Lo(NUV) = 3.9 x 10~ 14 x L (Ha), 
io(NUV) = 1.3 x Lo(FUV). 



(12) 



We adopt the factor k* = 1.4 x 10~ 28 x 7 in Eq. [T]for 
the NUV range and derive the values of modified fac- 
tors k* = 5.3 x 10~ 42 x 7 for the Ha emission line and 
k* = 1.8 X 10~ 28 x 7 for the FUV range. These values 
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correspond to the rough equality of the SFRo obtained 
from the initial Ha, FUV and NUV luminosities. The 
multiplier 7 can be used for overall tuning of the set 
of factors. It is equal to ~ 1, if the modified coeffi- 
cient fc*(NUV) i s set to its nominal value fc(NUV) by 
Kennicutt ( 19981) . On the othe r hand, if fc*(Ha) is set 
to its nominal value fc(Ha) by Kennicutt] ( 19981 ). then 
7 ~ 1.5. Thus, we obtain estimations of SFRo 

SFR = 5.4 x 1CT 42 x 7 x L(Ra)/f Ha (T) 

= 1.8 x 1CT 28 x 7 x L(F\JV)/f FUV (T) (13) 
= 1.4 x 10- 28 x 7 x L(NVV)/f NUV (T), 

where L(Ha), L(UV) and SFR are measured in erg s" 1 , 
erg s _1 Hz -1 and Mq yr -1 , respectively. The relations 
in Eq. [T3] give the approximately equal SFR (Ha), 
SFR (FUV) and SFR (NUV) for LCGs from our sam- 
ple. If the alternative correction for extinction is used, 
we obtain very similar results. Only one coefficient 
changes and the first row in Eq. [13] gets the form 
SFRo = 5.3 x 1CT 42 x 7 x L(Ra)/ f Ha (T). 

What value of the parameter 7 is preferable? We 
can measure the flux from galaxies in some wavelength 
ranges but we cannot directly measure their SFRs. 
These values can be estimated using different indica- 
tors of star formation. The values of the coefficients k 
in Eq. Q]for different wavelength ranges are obtained by 
modeling and therefore dependant from the parameters 
and assumptions in these models. Thus they can vary 
in some intervals. We use the certain values of k for the 
Ha line and the UV range. The values SFRo (Ha) and 
SFRo (FUV) are in general proportional, but not equal, 
as they must be. Thus we have to use some agreed set 
of the coefficients k for different ranges which matches 
various estimations of SFRo . We obtain the agreement 
condition for our sample in Eq. 1121 There is inter- 
val of the parameter 7 in which all coefficients k agree 
with values obtained from modeling. The increase of 
the quality of modeling will lead to fine-tuning of the 
value of 7. Unfortunately if we use the values of the 
coefficients in Eq. [T] for different wavelength ranges we 
get some discrepancy in SFRs obtained. This means 
that the models used to obtain these values have to be 
improved. In particular it would be useful to search out 
the ratio of initial values of the luminosities in differ- 
ent wavelength ranges to SFR after the onset of star 
formation. 



Lee et al.l fj2009h studied the consistency between 



the SFRs derived from the FUV continuum and Ha 
emission for a sample of the dwarf star-forming galax- 
ies. Particularly, authors discuss the dependence of the 
number ratio of ionising to non-ionising photons in the 
radiation of dwarf galaxies on its metallicity. We in- 
vestigated the dependence of both L(Ha)/L(FUV) and 



Lo(Ha)/Lo(FUV) on [O] and do not find any statisti- 
cally significant trend. All values of Fisher coefficients 
do not exceed F — 3.5. However, we note that the 
rather small range of [O] for the galaxies from our sam- 
ple makes it difficult to study this dependence. 

Using Eq. [T]we derive SFRs and obtain the distri- 
butions of SFR(Ha), SFR(FUV), SFR(NUV) as well 
as SRFo(Ha), SFR (FUV) and SFR (NUV). SFRs de- 
rived from the luminosities in the Ha emission line, 
FUV and NUV continuum vary in the wide ranges 
0.8 77 M© yr" 1 , 0.18 ~ 86 M Q yr" 1 and 0.24 113 
Mq yr -1 , respectively. The corresponding median val- 
ues of SFRs are 6.7 M yr" 1 , 3.8 M Q yr" 1 and 5.2 M Q 
y^ 1 . The median values of SFR (Ha), SFR (FUV) 
and SFRo (NUV) are 8.7 Mq yr" 1 , 5.1 M yr" 1 and 
6.5 Mq yr" 1 , respectively. For comparison, the me- 
dian value of SFR(Ha) is 0.92 M Q yr -1 for a sam- 
ple of about 7000 star-forming galaxies from the SDSS 
DR4 (galaxies being less lu m inous in H/3 than LCGs) 
( Izotova &: Parnovskv 20061 ). Cardamone et al. ( 2009 ) 
derived mean SFR ~10 M Q yr" 1 for the sample of 
"green peas" . 

Star formation rates SFRs derived from the Ha, 
FUV and NUV luminosities are in better mutual agree- 
ment if Eq. [13] is used instead of Eq. [1] We 
derive the median values of 5.8 x 7 x M yr" 1 for 
SFR (Ha), 6.6 x 7 x M yr" 1 for SFR (FUV), and 
6.5 x 7 x M© yr" 1 for SFR (NUV). 

We already mentioned in Section [4] that in general 
Eq. [T] for SFR can be applied for the continuous star 
formation during certain time interval At. Formally, 
for strongly star-bursting galaxies, we may estimate At 
as well, assuming the continuous star formation with 
the constant SFR and adopting nominal values for the 
coefficient k. Then, it is expressed as At = I/C5 with 
C5 from Table \5\ and attains the values in the range ~ 
5.1 - 8.9 Myr for different samples of case (b). 

The histogram for the Ha emission line luminosity 
L(Ha) is shown in Figure [71 The standard way to 
study the luminosity f unction is its ap proximation by 
the Schechter function dSchechterlll976T ) in the form 



t[j(L)dL = const(L/L*) a exp(-L/L*)d{L/L*) (14) 

where ip is the number of galaxies per unit volume in 
the luminosity interval from L to L + dL. Assuming 
that the volume V of the galaxies with luminosity L 
entering the sample V cx L 3 / 2 and using the maximum 
likelihood method we obtain the values a = —1.04 and 
L* = 8.5 x 10 41 ergs" 1 . This value of L* corresponds 
to SFR* = 6.7M yr" 1 according to Eq. [Q The value 
a = —1.04 is in agreement with the value obtained 
from 147,986 galaxy redshifts and fluxes from the SDSS 
( Blanton et~alll2003l ). 
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Fig. 7 Distribution of luminosities L(Ha) and SFR(Ha). 
Diagonal crosses show the distribution calculated from the 
Schechter luminosity function with parameters obtained 
from the maximum likelihood method. Error bars show the 
Poisson errors 

The distribution of galaxies calculated with the 
Schechter function is also shown in Figure [7l One can 
see that there is some difference between real and mod- 
elled distributions of Ha luminosities L(Ha). In partic- 
ular, model underpredicts the number of galaxies with 
largest luminosities calculated from Eq. [TJ] The dis- 
tribution of the observed luminosities in the range of 
low values near 10 41 erg s _1 is somewhat above the 
modelled one (see an inset in Figure [7]) . According to 
the Pearson's x 2 test this difference has the statistical 
significance over 99%. 

The distribution of Lo(Ha) differs from the Schechter 
function in a larger extent compared to the distribution 
of L(Hev) because of the larger luminosities. The distri- 
butions of the FUV and NUV luminosities are similar 
to that in Figure [7] after the luminosity scaling in ac- 
cordance with the difference of the median values. 



6 Summary 

We analyse the properties of Hcv and UV radiation 
for the sample of about 800 luminou s comp act galax- 
ies (LCGs) selected by llzotov et al. ( 201lh from the 
Data Release 7 (DR7) of the Sloan Digital Sky Sur- 
vey (SDSS). These low-metallicity galaxies are char- 
acterised by an active star formation with the star 



formation rate SFR(Hcv) in the range ~ 1 - 80 Mq 
yr _1 and can be considered as local counterparts of the 
high-redshift (z > 2) star-forming Lyman-break galax- 
ies (LBGs) and Ly-cv emission-line galaxies (LAEs). 
We use the optical SDSS spectroscopic data for LCGs 
to derive the luminosity in the Hcv emission line and 
SFR(Ha). Galaxy Evolution Explorer (GALEX) UV 
fluxes are used for obtaining luminosities and SFRs in 
the far-UV (FUV, A off = 1528A) and in the near-UV 
(NUV, A cff = 2271A) ranges. These data are supple- 
mented by o t her gl obal LCG characteristics derived by 
llzotov et al. ( 201 lh from their SDSS spectra: chemical 
element abundances of the interstellar medium, masses 
m and ages T of young starbursts. Our main results 
are as follows: 

1. We study the extinction in a sample of LCGs. 
It is found that LCGs are rather unobscured galaxies 
with a mean reddening about of E(B — V) = 0.136. 
The mean internal reddening in the sample is 0.106. 
For subsamples of "regular" galaxies with round shape 
and "irrregular" ones having shape with some sign of 
disturbed morphology we obtain internal reddening of 
0.081 and 0.120. The difference of internal reddening 
for subsamples is statistically significant value. 

2. We find that the ratio L(Hcv)/ra in starbursts 
with ages T < 3.2 Myr is constant implying that Hcv 
luminosity in young starbursts is proportional to the 
mass of the young stellar population. At later star- 
burst ages T > 3.2 Myr, the ratio L(Hcv)/m is declined 
exponentially with T. This temporal dependence of the 
L(Ha)/m ratio is in general agreem ent with that from 
the po pulation synthesis models by ISchaerer fc Vacca 
(1993) which predict the decrease of L(Hcv) after ~ 3 
Myr, the lifetime of the most massive stars. 

The dependences of the luminosities per unit mass 
of the young stellar population L(Hev)/m, L(FUV)/m 
and L(NUV)/m on T (Eqs. \5[ [9]) are similar implying 
that Hcv, FUV and NUV radiation is produced by the 
same young populations. However, the dependences of 
L(FUV)/m and L(NUV)/m on T are weaker as com- 
pared to L(Hcv)/to. Starting from To ~ 3.2 Myr, the 
half-times of the Hcv, FUV and NUV luminosities de- 
cline are 1.1 Myr, 1.6 Myr and 2.1 Myr, respectively. 
The ratios L(Hcv)/L(FUV) and L(Hcv)/L(NUV) also 
start to decrease after ~ 3.2 Myr. Thus, these ratios 
can be used for estimation of the starburst age T . 

With this value of T we can estimate m from L(Hcv) 
in young starbursts without invoking modelling of spec- 
tral energy distribution (SED). For that we introduce a 
function fn a (T) which takes into account the variation 
of L{Ha)/m with T. Then m - L(Ra)/ f Ha (T). 

3. The main impact of galaxy metallicity on its lu- 
minosity is the indirect one through the variation of 
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the mass of the young stellar population m. Dividing 
the sample of galaxies into subsamples with high and 
low metallicities, we obtain that the mean luminosity 
will be greater for the subsample with high metallicity 
due to the considerable increase of the mean value of 
m. On the other hand, the direct impact of metallicity 
has the opposite sign. The ratio L/m slightly decreases 
with increasing of the galaxy metallicity if the starburst 
age is constant. Thus, the galaxy luminosity decreases 
with increasing metallicity at fixed values of m and the 
starburst age T. This direct impact is weaker than the 
indirect one, but it is statistically significant. 

4. Luminosities in Ha and UV decrease rapidly af- 
ter the starburst age of ~ 3.2 Myr. We take into 
account this temporal evolution and introduce time- 
independent characteristics of the star formation activ- 
ity, namely the initial luminosities Lq at the starburst 
age T = 0. The initial luminosities in the Ha emis- 
sion line, FUV and NUV ranges can be obtained from 
the current luminosities and the starburst age T from 
Eq. HU We find that L (FUV) and L (NUV) are pro- 
portional to Lo(Ha) over a large range of luminosities. 
We can obtain the approximative equality of the val- 
ues SFR derived from the initial Ha, FUV and NUV 
luminosities by tuning the factor k in Eq [T] The set 
of factors k for Ha emission line and FUV and NUV 
ranges providing such equality for the sample of LCGs 
is used in Eq. [T31 

5. We find that SFRs derived from the extinction- 
corrected Ha, FUV and NUV luminosities vary in the 
wide ranges of 0.8^77 Mq yr" 1 , 0.18^86 M yr" 1 
and 0.24-^113 Mq yr" 1 , respectively. The correspond- 
ing median values of SFRs are 6.7 Mq yr -1 , 3.8 Mq 
yr _1 and 5.2 Mq yr . The median values of initial 
SFRs are SFR (Ha)=8.7 M yr"\ SFR (FUV)=5.1 
M Q yr" 1 and SFRo(NUV)=6.5 M Q yr" 1 . In all cases 
th e nominal coeffici ents fc(Ha), fc(FUV) and fc(NUV) 
by iKennicutt dl998h are adopted. The corresponding 
equalised median SFRo values in accordance with Eq. 
[Hare equal to 5.8 x 7 x M e yr" 1 for SFR (Ha), 6.6 x 
7 x Mq yr" 1 for SFR (FUV), and 6.5 x 7 x Mq yr" 1 
for SFR (NUV). 
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